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Abstract 

In this paper, a time substitution as used by Duru and Kleinert in their treatment of the hydrogen 
atom with path integrals is performed to price timer options under stochastic volatility models. 
We present general pricing formulas for both the perpetual timer call options and the finite time- 
horizon timer call options. These general results allow us to find closed-form pricing formulas for 
both the perpetual and the finite time- horizon timer options under the 3/2 stochastic volatility 
model as well as under the Heston stochastic volatility model. For the treatment of timer option 
under the 3/2 model we will rely on the path integral for the Morse potential, with the Heston 
model we will rely on the Kratzer potential. 
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I. INTRODUCTION 



Timer options, first introduced for sale by Societe Generale Corporate and Investment 
Banking (SG GIB) in 2007 jl, 2|, are relatively new products in the equity volatility market. 
The basic principle of this option is similar to the European vanilla option, with the key 
distinction being the uncertain expiration date. Rather than a fixed maturity time that is 
set at inception for the vanilla option, the expiry date of the timer option is a stopping 
time equal to the time needed for the realized variance of the underlying asset to reach a 
pre-specified level. 

Stopping times, sometimes formulated as first passage or hitting times, have applications 
in various research fields. Traditional applications of stopping times in physics are for 
example the situation where internal fluctuations induce the current of an electric circuit to 
a enseal value g Q a„d K.a.e.. p.oble. fl fl. Recent appl.caaon. of stop,,, 
times can be found in neuroscience, where a neuron emits a signal when its membrane voltage 
exceeds a certain threshold [t-O]; in the research field of quantum hitting times of Markov 
chains and hitting times of quantum random walks 
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111 ]: and in econophysics [l2|, Il3 |. 



or an introduction to first passage problems and an overview of possible apphcations see 



ll5|. 



When the expiration date is only determined by a stopping time that can theoretically 
become infinite the option is called a perpetual timer option. According to Hawkins and 
Krol [Q], it is usual practice to specify a maximum expiry for the timer option, at which 
point the option expires in the same manner as vanilla options, to prevent excessively long 
maturity times. These options are called finite time-horizon timer options. 

Timer options were first proposed in literature by Neuberger [l^] as "mileage" options in 
1990. In the middle 1990s, Bick emphasized the application of dynamic trading strategies 



with timer options to portfolio insurance as well as to hedging strategies [18| . Recently, after 
timer options were traded in the market, the amount of research concerning the pricing of 
perpetual timer options has increased. Li studied the pricing and hedging under the Heston 
stochastic volatility model jl9|. Bernard and Gui proposed a fast and accurate almost- 
exact simulation method in general stochastic volatility models 20|. Saunders developed 



an asymptotic approximation under fast mean-reverting stochastic volatility models 



We contribute to the existing literature by presenting analytical pricing results for both 
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perpetual and finite time-horizon timer options for a general stochastic process. These 
general results are then applied to determine explicit closed- form formulas for the 3/2 and the 
Heston stochastic volatility model. Especially for timer options it is relevant to investigate 
different stochastic volatility models, since the price of these options is particularly sensitive 
to the behavior of the volatility. 

We will derive these results in the path integral framework. The applications of path in- 
tegral methodology drawn from quantum mechanics to finance can be found, among others, 
3J27I and references therein. The virtue of this method lies in its capacity to provide 



m 



an intuitive way to derive the transition probability density function (propagator) of the 
underlying stochastic processes. Since the pricing of financial options comes down to eval- 
uating expectation values of stochastic processes, the path integral formalism is a suitable 
alternative to partial differential equations. 

To derive our general results we will rely on the Duru-Kleinert space-time substitution 



method used by Duru and Kleinert to treat the hydrogen atom with path integrals 22 
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291 ]. This method has recently been used in finance by Decamps and De Schepper to derive 



asymptotic formulas for Black-Scholes implied volatilities 30|]. The Duru-Kleinert space- 
time substitution approach serves here to translate the original stochastic processes to new 
ones behaving in a stochastic time horizon. Under this new time horizon, the random expiry 
time is expressed as a functional of the transformed stochastic volatility. Then a method 
related to variational perturbation theory 22l, |3l| is applied to derive the joint propagator of 
the transformed stochastic volatility process and the stopping time process in the new time 
horizon. Based on these transition probability density functions, we arrive at the pricing 
formulas for the perpetual timer option. In addition, we obtain pricing formulas for the 
finite time-horizon timer option by deriving the joint propagator of the log-return and the 
realized variance process. 

For stochastic volatility processes, we start by emphasizing the 3/2 stochastic volatility 
model 32-3^ not only in view of its analytical tractability but also because of the support 
from empirical evidence 35l-l37j]. The results for this model are obtained by making a 
connection with the Morse potential. Next we treat the Heston stochastic volatility model 
by relating it to the Kratzer potential. This leads to closed-form pricing formulas for 



perpetual and finite time-horizon timer options for both models. The result for the perpetual 
timer option under the Heston model corresponds to the one found by Li [l9^, confirming 
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our approach. 

This paper is organized as follows. In section [TTl we present general pricing formulas 
for perpetual and finite time-horizon timer call options under general stochastic volatility 
models. Section UTTl is devoted to deriving closed-form formulas for the 3/2 and the Hes- 
ton stochastic volatility model. In section IIVI the closed-form formulas are compared with 
Monte Carlo simulations and some properties of timer options are discussed. And finally a 
conclusion is given in section |Vl 



II. GENERAL PRICING FORMULA OF TIMER CALL OPTIONS 



A. Model description 

For conciseness of representation, in this paper we only consider option pricing in a risk- 
neutral world. Moreover we assume that the initial time of the option is the current time 
t = because the generalization to the case of a forward-start option is straightforward. 

Let {5'(t)} denote the underlying asset price process following a Black-Scholes type 
stochastic differential equation (SDE), with a variance v{t), which is stochastic variable 
itself. Conventionally, the time evolution of S{t) is represented in terms of the log-return 
x{t) = In^^, with 5*0 = S{0). The realization of a stochastic process Z at a special time 
s will be denoted by Zg, and we will use this notation throughout the paper. After the 
transformation to the log-return, the system is governed by the SDEs: 

dx{t) = (^r - ^) rft + (Vl - pWi + pdW^ , (1) 
dv{t) = a{v)dt + l3{v)dW2, (2) 

where r is the constant risk- neutral interest rate, Wi{t) and W2{t) are two independent 
Wiener processes, p G [—1, 1] is the correlation coefficient between x{t) and its variance v{t). 
Now we introduce the notion of the realized variance, which is a principal ingredient of 

N 2 

timer options. In practice the realized variance is given by ^ {xt^ — Xt„_j^) , where the set 



n=l 

of evaluation times are for example daily closing times. In the literature (see 19|, |39| ) . 



the realized variance of the underlying asset during a time period [0,T], denoted by It, is 
usually approximated by: 

It = [ v{t)dt. (3) 



Also in this paper equation ([3]) will be used as the definition of the realized variance. 



B. Pricing of perpetual timer options 

The price of a perpetual timer call option with strike price K can be expressed as the 
expectation of the discounted payoff: 

Cperp = E [e-'^" max {So e^^e - K, 0)] . (4) 

This expression is similar to the one for the vanilla call option, except for the uncertain 
expiry time Tg, which is the stopping time defined as 



|m > 0; ^ v{t)dt = . (5) 



Tb = inf |m > 0; J v{t)dt = B 

Here B = ct^Tq is the pre-specified variance budget with Tq the expected investment horizon 
and (To the forecasted volatility of the underlying asset during that period. 

The dependence on the implicitly defined expiry time Tb is inconvenient. We will now 



apply the Duru-Kleinert method of quantum mechanics 22| to construct variables in function 



of which Tb is explicitly given. Define a time substitution r(t) such that 



r(t) = / v{s)ds, (6) 







we will refer to r as the pseudotime, following 
that 



22| . The inverse function theorem gives us 

dT-\t) 1 



dt v{T-^{t))' 
from which it follows that T^^(t) is given by: 



(7) 
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Denote f (r~^(t)) by V{t) and x{T~^{t)) by X{t), which follow new SDEs: 

dVit)=^dt + ^dW,, (9) 
dX{t) = (^^ - rft + - p'dWr + pdW^ . (10) 

Given the timer call variance budget 

B= r v{t)dt = T{TB). (11) 

Jo 



we obtain the explicit expression for the stopping time as 



Tb 



T-\B) = / 
Jo 



-dt 



Vit) 



dt. 



(12) 



Note that {x{t),v{t)) = (X(r(t)), F(r(t))) , so as {x{t),v{t)) evolves in the period [0,7^], 
{X{t),V(t)) evolves in [r(0), r(7B)], that is [0,B]. Therefore [0,B] is now a fixed horizon in 
pseudotime, and not only do the processes X, V and T evolve during that period, but also 
expression (j4j) can be written as 



Cperp = E [e-^'^" max (^o - K,0)] . 



(13) 



Hence, it is intuitive to study the joint transition probability density function of the dynamics 
of (X, T). However, as T depends on V, we turn to the joint propagator of the dynamics 

of (x,\/,r). 

The substitutions 

w 



z{t) 



y{t) =X{t)-pz{t) 



dV{t), 



(14) 
(15) 



help change the correlated dynamics of {X, V) into two independent processes following 



dy{t) 



[V{z) ^ 
dz{t) =A (z) dt + dW2 



pA{z) 



dt + v^l - p^dWu 



where 



A{z{t)) 



a{V) 1 d f ^A^\ P'^{V) 



+ t: 



V 



(16) 
(17) 

(18) 



p{v)Vv ■ 2dV \PiV) J 
is a function of z(t) because V{t) is expressed in terms of z(t) according to expression ([1] 

To determine the price of the timer option, the propagator V {ye, zbjTbIvo, zo,0) is 
needed. This propagator describes the joint probability that y has the value yts, z has 
the value zb and the stopping time has the value 7b at a later pseudotime B given their 
initial value yo, zq and at pseudotime 0. Since the processes y and z are uncorrelated, the 
Lagrangian corresponding to their joint evolution can be written as £[?/, z\ + £[2;, i] with: 



C[z,z\ 



y 



V(z) 



pA {z) 



2 1 



^Wf + i|^(.), 



(19) 
(20) 
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Using the path integral framework, the joint propagator V {i/b, zb, Tb \ Ho, zq, 0) can be de- 
termined by: 

V {vb, Zb, Tb \ yo, zq, 0) = 

JvyJ Vz6 (Tb - J^^ v^^^) e-/o"(^[^'^l+^b.^.-l)'^*, (2i) 

where 6{-) is the deha function. It serves here to select these paths of V, expressed in terms 
of z, such that ^^dt equals Tb- 

To proceed we introduce the Fourier transform of the delta function. Furthermore since 
the path integral corresponding to the y variable is quadratic it can be solved analytically. 
After performing this path integral, we can return to the original Xb variable. Expression 
fl2T|) then becomes: 

V {XB,ZB,TB\yo,zo,0) = 

_[Xs-T(z)]2 

,27r J v^27r(l-p2)i3' ^ ^ 



where 



T{z)=p(^ZB-zo- A {z{t)) dt^ + tTb - |. (23) 

In order to add the z dependent term T {z) to the Lagrangian of the z path integral one can 
introduce another Fourier integral, and V (Xg, zg, Tb \ yo, zq, 0) then becomes: 

V {Xb, Zb, Tb \ yo, zo, 0) 

oo t; (i-o'^']b poo 



dl^^lx,-^^l^ r dp 
27r 2n 



X [ Vz{t)e'-^°^^^"'"^^'''^^''^^^"^^'^\ (24) 



Whether our approach will lead to closed-form pricing formulas for timer options will depend 
on the Lagrangian C[z, z] + ipyj^ + H'^ {z). More precisely this means that A {z{t)) and 
T {z) should be well behaved enough in terms of z. For the two examples illustrated in this 
paper, the functions a{y{t)) and (5{V{t)) are as such that T {z) is only a function of zb 
and Tb, denoted by T {zb, Tb)- Then it is not necessary to introduce the Fourier transform 
of expression (1211) and we can proceed with expression (12^ . Now the price of a perpetual 
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timer option which is given by 

Cperp = j dXs j dZB j dTs V (Xb, Zb, Tb I 1/0, Zo, 0) 

X [e-'"^« max (Sq - i^, O)] , (25) 

can be written as: 

/•OO /"OO 

Cperp = dTB dzBV {zb, Tb I zo, 0) C {zb, Tb) , (26) 

Jo J-oo 

with C {zb, Tb) being the prices conditional on z: 
C{zb,Tb) 

[Xs-r{zBTs)f 

/°0 2(l-p^)S 
rfX^ [e-^^« max (5o - K, O)] 

-OO A/27r (1 - p^) i3 

=^^eT(.H,r8)-rr8+(i^_^(^^) - K e~'^''Afid^), (27) 



where A/'(-) is the cumulative distribution for the normal random variable and 

_ln% + {l-p^)B + TizB,TB) 



Inf + T(^6,rB) 



(28) 
(29) 



and P (zz3, Tb \ zq, 0) is given by 

V{zb,Tb\zo,0) = r ^e'^^^ [ p^e-^'^(^t^''l+^)'*. (30) 

Note that expression (127|) is a Black-Scholes-Merton type pricing formula for perpetual 
timer options. To determine the price of a perpetual timer option for a particular model one 
needs to evaluate T {zb,Tb) in order to obtain C {zb,Tb)- Furthermore if we also have the 
analytical expression for the joint propagator V {zb,Tb \ Zo,0), formula (126|) demonstrates 
that the closed-form perpetual timer option pricing formula can be derived through two 
trivial integrals. 



C. Pricing of finite time- horizon timer options 

In this subsection we consider the pricing of finite time-horizon timer option. 
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Let the maximum expiry time to be T, then the price of a finite time-horizon timer 
option, denoted by Cpini, of strike price K can be expressed as a sum of two contributions: 

Cpini = Ci + C2, (31) 

where 

Ci= dTB dzsVizts,%\zo,0)Cizs,rB) (32) 

Jo J -00 

is the contribution from paths that exhausted their variance budget before time T, and 

/oo 
(^0 e"^ - K)^ Vb (xt \xo ) dxr (33) 
-00 

is the contribution from paths that reach the preset finite time horizon. Note the integration 
range of Tb in Ci, which is truncated by the maximum expiry time T. 

Denote the joint propagator of the log-return and the reahzed variance as V {xt, It I Xq, 0), 
then 

Vb {xt \xo) = V (xt, It I Xq, 0) dlr- (34) 
Jo 

Note Vb {xt \ Xq) is not the propagator of x that should be used for the European vanilla 
option, which represents the probability that x has the value xt at later time T given the 
initial values xq at time 0. Instead Vb {xt \ xq) in C2 is the propagator of x which is also 
conditioned on the fact that the realized variance budget of each path has not been exhausted 
before the maximum expiry time T. 

Furthermore, if Vb {xt \ xo) can be written as a Fourier integral: 

^^Ui^T-rT)j,^^^^ (35) 



then by following the derivation outlined in 40|, we can rewrite C2 explicitly, and thus the 



pricing formula of finite time-horizon Timer option as 

Cpini = dTB dzBV {zb, Tb \ zq, 0) C {zb, Tb) 



JO J -00 

5(0) ^ . r dl e"i"'t-.-'^)g{i.) 

where 

g{l) = S^T{l + i)-Ke-''''T{l). (37) 

The integration of Vb{xt\xq) over all possible xt's gives v^(0), which is the "sur- 
vival probability" describing the probability that the underling asset is executed at 



the maximum expiry time T. This survival probabihty can also be determined by 
dTs J dzjsV {zs,Tb\ zo,0), from which it is clear that this probability is independent 
from the evolution of x, thus does not depend on the correlation coefficient p. 

For finite time-horizon timer option, besides the evaluation of propagator V [zjs, Tb \ Zq, 0) 
as for the perpetual timer option, we must also calculate the propagator V (xy, It\xo, 0) to 
derive the formula of J-'{1)- 



III. PROPAGATORS FOR THE 3/2 AND THE HESTON MODEL 

In this section we focus on the derivations of joint propagators V (zg, Tb \ zq, 0) and 
V {xt, It I xq, 0). These are used in section [IV] in conjunction with expressions (!26|) and (!36|) 
from the previous section to price perpetual and finite time-horizon timer options, respec- 
tively. Note that V {zb, Tb \ zq, 0) is evaluated in the pseudotime horizon and V {xt, It \ xq, 0) 
in the original time horizon. The 3/2 and the Heston model are chosen both from mathe- 
matical and empirical considerations. 

As mentioned in the previous section, it is convenient to choose models such that 
Jq a {z{t)) dt is a function of zb and Tb- In addition, from the perspective of mathematics, 
the total Lagrangian in expression (!30|) : 



CTot[z,z]=C[z,z] + ^ (38) 

written in terms of z and z should be well behaved enough to achieve a closed-form solution 
with the path integral. 

Furthermore there is substantial empirical evidence supporting the stochastic differential 
equation underlying the 3/2 model. The Heston model, on the other hand, is important 
because it is a standard model for the financial industry. 

A. The 3/2 model and the Morse potential 



The model dynamics of the 3/2 stochastic volatility model [32| is given by: 



dv{t) = Kv [6 — V ) dt + ev^^^dW2. (39) 
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Relating this model to the general stochastic volatility model used in (j2]), we have 

a iV) =kV (e-V), (40) 
/3 (V) =eV^/\ (41) 

For calculation convenience, we multiply z{t) defined in expression f lT4|) by a factor — e to 



obtain 



z{t) = -\nV{t). 



Thus, according to equations ffT^ . f lTS]) and f l2^ . we have 



B 



A{z{t))dt 
T (zb, Tb) 





—Tb 
e 



[zb + In Vq) + rTB- — 



(42) 

(43) 
(44) 

(45) 



Therefore the total Lagrangian is 
^Tot [z 



1 .2 






+ ^e'z 


-(5 



,22 



+ 2 



— ^ + — «p ) e 



2e2 



(46) 



The nontrivial terms of Cxotaiz^z] reveal that z{t) is subjected to a Morse potential. By 
making use of the known path integral for the Morse potential see Appendix |X1 the 
joint propagator is expressed as 

V {zb, % I zq, 0) 



dp 
2^ 



,2\2 



e /"oo 



dp 
2^ 



X e^^'7-B|p^(^)^-/o1^ + #e--(^+«.-.p)e^].* 
'^^ -ii^(e^S-e-0)+(^ + i)(^g_2o) 



e2 sinh ^ 



e £" 



X e '^I^'^aJ 2 

d<l> 



B+ ( ^ + 1 ) KeTs - (e^B +e^0 ) coth 



X 







-Re 



;2- C 2 



V^* I sinh 



(47) 
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with J (■) the modified Bessel function of the first kind, and under the condition that the 
real part of the integration variable satisfies 

2 



2 1 



(48) 



according to (IA6p . 

Plugging expressions ( H5ll and ( H71) into formula ( l26l) yields the closed-form pricing formula 
for the perpetual timer call options under the 3/2 model. 

The integral over all possible can be done analytically, which leads to the marginal 
propagator for the stopping time 7^: 



K0 



kOTb 



V {% I 0) =^ ( 1 + coth \ e-(^+l)'4^ 



XI — ^ e^" I — 

'0 



TT 



r 2,/4$ + i 



X ,Fl|>^ + l;2^/4$ + l;-^l, 



(49) 



where F (■) is the Euler gamma function, iFi (■; ■; ■) is the confiuent hypergeometric function, 
and 



J^iTs) = ^ (coth 



(50) 
(51) 

We move on to the calculation of the propagator V {xt, It I ^^o, 0) by performing the 
following substitutions 



x{t) = X — - {Inv — K6t) — rt, 
e 

m = ^, 



(52) 
(53) 



which lead to two uncorrelated processes 
dx{t) = 

dm = 



k6 ^ f K 3 n \ 1 



dt - -dW2, 
2 



(54) 
(55) 



12 



P(V,.T„|V„,0) 





0.05 0.1 0.15 0.2 



FIG. 1: (Color online, two columns) This figure shows several aspects of the time evolution of 
variables relevant for timer options under the 3/2 stochastic volatility model. The upper left panel 
shows several simulated variance paths up to the point where the realized variance reached B. The 
bottom left panel shows corresponding log-return paths. The inset of the upper left panel shows 
the probability distribution of the stopping time Tb- The inset of the bottom left panel shows 
the density VBi^Tlxo) determined by expression The upper right panel shows the joint 

probability distribution of the variance and the stopping time. The bottom right panel shows the 
joint probability distribution of the log-return and the realized variance. The parameters used here 
are: vq = (0.295)2, k = 22.84,6* = (0.4669)2, e = 8.56, S = vo,r = 0.015,/) = -0.5, T = 1.5. 



and thus the corresponding Lagrangians: 



2v (1 - p2) 
£[(,(] = C,[C] + C,[Q, 



(56) 
(57) 
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where 



A[C,C] 



4 



+ 



2Ke- 



2k 3 



c2 



(58) 



(59) 



2 

Since x is independent from ^, so the probabihty that x go^s to xt, C goes to Ct and the 
reahzed variance reaches It at a later time T given the original positions xoi Co and /q = 
at the initial time is 



^(Xt, Ct, I Xo, Co, 0) 



T 



VC{t) ^yr- j v{t) dt ] e-/o ^[^'^H* 



2k I 3 



C( 



+00 



dl 

2^ 



]^Ct 
Co 



— oo 



dp 
2^ 



X / Vat)e 



Jo 



2m2 



.Ci[C,C]+- 



(60) 



where the remaining path integral over T>({t) of the radial harmonic oscillator potential 41 1 
given by: 

2Key/CTCo ^-«|fc?.+Cf^)coth(^) J / 2k9 CtCo 



with 



sinh(- 2 j 



sinh(- 



«2 



2/ (-i + f + f)+^i^ 



+ ip 



(61) 



(62) 



Integrating over Ct leads to V (xt, It \ Xo, 0): 

^ (Xt, I Xo, 0) 

/■oo 

= / 'P(Xt, Ct, I Xo, Co, 0) dCr 
Jo 



dp 
2^ 



+ 00 77 / 9 



2tx 



M 



e2Ar 
2 

'1^ 



(63) 
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where 



(64) 
(65) 



Expression fl63l) agrees with expression (73) in [37 1 
According to (134|) . we have for the 3/2 model 



/+00 77 
« e./(x.-rT)^^^^^ (66) 
-oo 27r 



where 



T(l\ -_•/"" e'^^ - 1 / r(A + 1 - M) 



r(A + i) 



X + (67) 

The closed-form pricing formula for the finite time-horizon timer call options is derived 
by substituting ( !66l) and ( !67l) in expression (!36l) . 

B. The Heston model and the Kratzer potential 



For Heston stochastic volatility model 38|], the model dynamics is written as: 



dv{t) = k[9 — V ) dt + ay/vdW2. (68) 

To relate this model to the general stochastic volatility model (jj]), a {V) and f3 {V) are given 
by 

a{V) = K{e-V), (69) 
I3{V) = aW. (70) 

From equation f|T4l) . we have the relation between z{t) and V{t): 

z{t) = ^V{t). (71) 

thus the stopping time is a functional of z(t): 

1 /-^ 1 



= - / —.dt. (72) 
Jo 4t) 
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Plugging equations (168|) and ([2]) into definition (fTSll gives 

r A{z{t))dt = —TB--B, 
Jo a a 

therefore (written in original variable fo) 



(73) 



T (zb, Tb) =P (zb Tb + -B 

1 — i 
£toJ2;, il =-z^ H - 

1\ i 



tTb - 



(74) 



where 



A 



(75) 



(76) 



The nontrivial terms of the total Lagrangian CTot[z-,z\ manifest that z{t) is subjected to a 
Kratzer potential. With the help of the known path integral for Kratzer potential [dJ], see 
Appendix [Bl we obtain the joint propagator as: 



V {zb, Tb I zo, 0) 

V ^io z{t) J 



dp it 



ipTe 



27r 



X / Vz{t)e 



Jo 



1 -2 I ~1 



(^+^)a'-^ 



2 4 



X Re 



TT 



g<I>B ^-v25(zB+2o) coth(y -2<^'7e 



X ■ 



2V25 r 

-'2A 



sinh( "2 <^Te 



sinlWy^crTe 



(77) 



According to the condition given by (IBSp . the real part of the complex variable $ must 
satisfy 



K 



2cr2 



(78) 



Again, plugging the expressions ( !7i|) and (!77|) into formula ( 126|) yields the closed-form pricing 
formula for the perpetual timer call options under the Heston model . 
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Note the similarity of expression (177|) with the result obtained in 19|]. Expression (177|) 
has a clear meaning as joint transition probability density function which illustrates the 
conciseness of physics; path integration allows to derive expression (1771) without applying 
any previous results of Bessel processes as done in [l^. 

Li [l9| computed the risk-neutral expected maturity in expression (5.2) by doing two 



numerical integrals. Actually, we can derive the marginal propagator of the stopping time 
Tb by integrating over all possible as follows 



V{Tb\^) 

poo 

I V {zi3,TB\Zo,0)dzB 

Jo 



'0 



IT 



X Re 



exp <^B + 



2A+1 



X 



( 



2A + 1 + 



2A+2 



/i+V2¥coth(A/f o-Te 



v 



\ 



(79) 



For the calculation of V {xt, -^t|0, 0), we follow the derivation in [25|. Substitutions 



x(t) = X (f — nOt) — rt, 



a 



(80) 
(81) 



give two uncorrelated processes: 



pK 1 



dxit) = - - ]vdt + v^Vl - P^dWu 



dm - 

The corresponding Lagrangians are 



2C 2^ 



1 



2v (1 - p2 

£[C,C] = A[C,C] + A[C,C] 



dt + -dW2. 



PK 1 



(82) 
(83) 



(84) 
(85) 
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where 



£-2 



2k9 1\ C 2k 



2J C ' ■ 



(86) 
(87) 



Since x is independent from we similarly have the joint propagator of the dynamics of 
X, C and J: 



'PiXT, Ct, It I Xo, Co, 0) 



vat) 5 [It - 



X / Vx{t) e- ^^X'XMdt 



2kB_1 

Co 



dp ■ 



X 



27T 



X / PC(t)e 



- r 

Jo 



dt 



(88) 



where the path integral for the radial harmonic oscillator potential 411] is given by: 

Au CtCo 



0"^ sinh(ci;T) 



2lj 



(C|.+Co')coth(a;T) 



•--\ 



0"^ sinh(a;T) 



where 



CT / 



+ (l-p2)/2+^/ 



2 y ^2 

Integrating over Ct leads to V (xt, -^t | Xo, 0): 



2pK 



a 



2ip. 



(89) 



(90) 



^ (Xt, It I Xo, 0) 

= / VixT, Ct, It \xo, Co, 0)dCT 
Jo 



.^^vo+^T 



dp 
2^ 



dl , 



U{xT-rT) 



2tt 



2kB _ 2i^(coBh(^T)-iV) , 



(91) 



where 



= ( cosh(a;T) H ^ sinh(ci;T) 



(92) 
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FIG. 2: (Color online, two columns) This figure is similar to figure 1, but now for the Heston 
stochastic volatility model. The parameters used here are: vq = 0.087, k = 2,9 = 0.09, a = 
0.375, e = vo,r = 0.015, p = -0.5, T = 1.5. 



Note the similarity of expression (l9Ti) with the result obtained in 34 1. 
According to fl34|) . we have for the Heston model 

/ + 00 TT 
^^Ui^T-rT)j,^l^^ (93) 
-oo 

where 

dp e'P'^ - 1 



271 p 



^ ^iiE(^eT+vo)j^^ g — .^sinh(^T) (94) 

with which we obtain the closed-form pricing formula for the finite time-horizon timer call 
options according to formula fl5B]) . 
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TABLE I: Comparison of the analytical and the Monte Carlo (MC) simulation values for both 
perpetual (Cperp) and finite time-horizon {Cpini) timer call option prices under the 3/2 model. 
The columns indicated by RE shows the relative error (in %). Parameters used here are: vq = 
(0.295)2, K = 22.84,61 = (0.4669)2, e = 8.56, = vo,r = 0.015, T = 1.5. 



K 


P 


^Perp 




Analytic MC RE(%) 


Analytic MC RE(%) 


90 


- 0.5 


0.5 


17.8064 17.8128 -0.0359 
17.7046 17.7129 -0.0469 
17.5839 17.5853 -0.0080 


17.6813 17.6790 0.0130 
17.5385 17.5510 -0.0712 
17.4260 17.4301 -0.0235 


100 


- 0.5 


0.5 


12.5780 12.5784 -0.0032 
12.4619 12.4683 -0.0513 
12.3300 12.3231 0.0560 


12.4089 12.3998 0.0734 
12.2780 12.2890 -0.0895 
12.2104 12.2032 0.0590 


110 


- 0.5 



0.5 


8.6518 8.6414 0.0486 
8.5339 8.5388 -0.0574 
8.4026 8.3943 0.0989 


8.4381 8.4301 0.0949 
8.3531 8.3611 -0.0957 
8.3229 8.3153 0.0914 



IV. PRICING RESULTS AND DISCUSSION 



In the previous sections explicit formulas concerning timer options are derived for the 
3/2 and the Heston model. The analytical tractability of these formulas is demonstrated by 
figures [1] and [2] and tables I and II. Figured] as well as table I are devoted to the 3/2 model 
and figure [2] as well as table II to the Heston model. For the 3/2 model the parameters are 



based on Ref. 



34j |. where they are calibrated on market prices of S&P500 European options 



to guarantee the relevance of these parameters. The parameters for the Heston model were 
chosen such that the two models are comparable. 

The two tables compare timer option prices calculated with the formulas of the previous 
sections with prices obtained by Monte Carlo simulations. Results for both the perpetual 
and the finite time horizon timer option are presented for several strikes and correlation 
values. For all the Monte Carlo simulations presented here, we used 20 million samples and 
3200 time steps per year. The relative error between the exact and the simulated prices is 
always less than 0.1% confirming our formulas. 

Although the prices of the timer option presented in these tables vary only slightly as 
a function of the correlation coefficient p, the timer option does have different features for 
different correlation values. This can be seen in figures [1] and [2l The upper left panel shows 
possible realizations of the variance up to the point where the realized variance reached B. 
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TABLE II: Comparison of the analytical and the MC simulation values for timer call option prices 
under the Heston model. The columns indicated by RE shows the relative error (in %). Parameters 
used here are: vq = 0.087, k = 2, 6* = 0.09, a = 0.375, B = vo,r = 0.015, p = -0.5, T = 1.5. 



K 


P 


^Perp 




Analytic MC RE(%) 


Analytic MC RE(7o) 


90 


- 0.5 


0.5 


17.8095 17.7948 0.0826 
17.7249 17.7232 0.0096 
17.6263 17.6146 0.0664 


17.6914 17.6851 0.0356 
17.5351 17.5330 0.0120 
17.4627 17.4680 -0.0303 


100 


- 0.5 


0.5 


12.5789 12.5668 0.0963 
12.4772 12.4763 0.0072 
12.3691 12.3586 0.0842 


12.4034 12.4010 0.0194 
12.2675 12.2678 -0.0024 
12.2426 12.2464 -0.0310 


110 


- 0.5 


0.5 


8.6515 8.6412 0.1192 
8.5449 8.5446 0.0035 
8.4393 8.4317 0.0890 


8.4206 8.4218 -0.0142 
8.3393 8.3405 -0.0144 
8.3522 8.3542 -0.0239 



The lower left panel shows the corresponding log-returns. In these figures we used a negative 
correlation. As a consequence, paths with a low log-return are more likely to have a high 
volatility and the corresponding option will probably be exercised sooner than an option 
with a high log-return. This behavior is also seen in the inset of the lower left panel. This 
inset shows the density given by formula ( IMl) when the maximum expiry time T reached 
1.5. Recall that this is the distribution of log-returns whose realized variance has not yet 
reached B. Due to the negative correlation, paths with a low log-return are more likely 
to have reached B and will therefore less likely contribute to this distribution than paths 
with a high log-return. Therefore this distribution is clearly shifted to the right. The lower 
right panel shows the joint density of the log-return xt and the realized variance It when 
T equals 1.5 and also illustrates this behavior. 

Tables I and II illustrate that the prices for timer options are quite similar for the two 
stochastic volatility models. Nevertheless there are important differences between the two 
models concerning timer options. This is illustrated by the upper panels of figures [1] and 
[21 As already mentioned, the upper left panel shows several possible time evolutions of the 
variance. The inset of this panel shows the probability distribution of the stopping time Tb- 
The upper right panel shows the joint density of the variance and the stopping time, which is 
useful for an intuitive understanding of the time evolution of the underlying processes. The 
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explicit form of this density is not included in the text because it is not needed to calculate 
prices and because it can easily be derived from expression fl30|) . For the 3/2 model, the 
probability that the variance reaches large values is larger than for the Heston model, while 
the probability that the variance reaches very small values is smaller than for the Heston 
model. Therefore the probability that the timer option will be exercised very fast is larger 
for the 3/2 model than for the Heston model. On the other hand for the Heston model there 
is a larger probability that the timer option will only be exercised after a long time than for 
the 3/2 model. 



V. CONCLUSION 



In this paper we construct a method to price both the perpetual and the finite time- 
horizon timer option for a general stochastic volatility model. Pricing of such options is 
related to first passage time problems in that the stopping time for the option is determined 
by a boundary on a cumulative stochastic process. The method proposed here is based 
on the Duru-Kleinert time transformation and the path integral framework. Furthermore 
we discuss the conditions a stochastic volatility model has to satisfy in order to be able 
to derive closed-form pricing formulas. These general results are then applied to derive 
closed- form formulas for the Heston and the 3/2 stochastic volatility model. For the 3/2 
model this involves the solution of the Morse potential, for the Heston model the Kratzer 
potential needs to be solved. Finally, our closed-form pricing formulas are shown to be 
computationally tractable and are validated by Monte Carlo simulation. 



Appendix A: Path integral for Morse potential 



The path integral for the Morse potential is given in Ref. 41| : 



2 2m V / 



dt 



mV (i + V-2mE/h - aVo) 

n?VoT (l + 2V-2mE) ""Vb,/ 



W 



2m E 



{2Voe 



„ T / / 2 m J 



i2Voe^< 



(Al) 
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where Vq > 0, x> (x<) is the larger (smaller) of two variables x' and x", r(-) is the gamma 
function, and M...(-) and are the Whittaker functions related to the Kummer and 

Tricomi confluent hypergeometric functions, respectively. 

This path integral (let h = 1, —E = $ = $^ + i<l>j), analytically continued to "imaginary 
time" defined by 

t = -ir (r G M), T = -iA (A G M), (A2) 

results in: 



dAe-"^ I Vx{T)e 



- 

JO 



2 \dT } "'"2m 



-{x'+x")/2 



VqV {l + 2V2m$j 

(2yoe->)M,^^_^ 



(2V^oe"^ 



(A3) 



Therefore the inverse Laplace transform gives the expression of the propagator evolving in 
the "imaginary" time horizon [0, A] as shown in Ref. |22|: 



Jo 



f(f)' + ^(-'^-2-^) 



dr 



r (^i + V2m$ - aVo 
VoT (^1 + 2V2m$j 

^^aVo,V^^ 



Moo ^^^^m 



e 2 



2m — ^Re 
Jo TT 



-Vb(e^'+e^") coth^ ^ 



sinh ^ 

2Vo e 2 



sinh^ 



under the condition 



that is 



Re 



+ V2m$ - aVo 



> 0, 



{Re[a]Vo-^y 



2m 



(A4) 



(A5) 



(A6) 
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Appendix B: Path integral for Kratzer potential 



The path integral for Kratzer potential is given in Ref. 4l| (A > 0, k = ^^^/~^)- 



■? / ™ / ' f*" / ™y.2 1 £l£2 r? ^ -f 1 J, 



h 



dTe'''^/'' / Vr{t)e 



1 ri^T{l + X-K) 



hy 2E r(2A + i) 

X W^^x (v/=8^ y ) M,,A (v^=8^ y ) , (Bl) 

The analytical continuation of this path integral (let h = l,m = 1, 6162 = /3, —E = $ = 
$ij + as used in Appendix El gives ( k = 



r{A)=r" 



c?Ae 



Vr(T)e 



r(0)=r' 



r(i + A-«:) 



/2<i>r(2A + i: ^"'^ ■ ^"'^ ■ ^' ^^^^ 

Therefore the inverse Laplace transform leads to the propagator in time horizon [0, A] as in 
Ref. |42l: 



Pr(r)e 
=° d<l>i 



Jo 



1 Mr\2 A^-l ^ 



Re 







TT 

2/35-V2$(r>+r<) coth(V2*?) 



X e 



yrvrv / 2V2<I> ^r>r< 

X 7 ^ J2A 



sinh(y2$n \sinh(y2$C 



valid under the condition 



that is 



Re 



2 



> 0, 



1 / /3 



2 V A + 



(B3) 

(B4) 
(B5) 
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